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Abstract. A nonlinear development of radial pulsation instability to a 
resonant excitation of nonradial modes is studied. The theory covers the 
cases of axisymmetric (m = 0) modes as well as (m, —m) pairs. Adopting 
a simplified treatment of the mode coupling, it is found that multimode 
pulsation with constant amplitudes is a highly probable solution. An 
observable consequence of the m = mode excitation is randomness of 
pulsation amplitude. The case of an I = 1 mode is the most important 
because of a small averaging effect and a large excitation probability. A 
significant amplitude and phase modulation is predicted in the case of 
excitation of an I = l,m = ±1 pair. This may explain Blazhko-type 
modulation in RR Lyrae stars. If this model is correct, the modulation 
period is determined by the rotation rate and the Brunt- Vaisala frequency 
in the deepest part of the radiative envelope. 



1. Introduction 

One of the most intriguing questions concerning RR Lyrae stars is the nature of 
the long-term modulations that characterize about 20-30% of RRab and 2-3% 
of RRc stars (e.g., Moskalik & Poretti, these proceedings; Kovacs, 2000). The 
modulations of the oscillation amplitudes and phases manifest themselves in the 
Fourier spectra as secondary peaks with frequencies close to the main pulsation 
frequency. There seems to exist two distinguishable cases: one secondary peak 
with frequency either higher or lower than the main frequency and a pair of side 
peaks which, together with the main peak, form an equidistant triplet. 

There are two models proposed to explain the modulations: a magnetic 
oblique rotator (Shibahashi, 2000) and the 1:1 resonance model (Van Hoolst et 
al., 1998). In this paper we develop the resonant model. 



2. Nonresonant properties of oscillation modes 

As was shown by Van Hoolst et al. (1998) and Dziembowski & Cassisi (1999), 
there always exists a very dense spectrum of \ow-£ nonradial modes in the fre- 
quency range of the lowest radial overtones in models of RR Lyrae stars. The 
most unstable modes have frequencies very close to the fundamental as well as 
the first overtone frequencies. Among the low degree modes in the closest vicin- 
ity of the radial modes, the I = 1 modes are the most unstable and thus the 
most important dynamically. 
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Small frequency differences between the radial and the closest nonradial 
mode suggest that 1:1 and similar resonances should be taken into account in 
the study of the oscillations. Moreover, this frequency difference, which is an 
important parameter in resonant dynamics, is such a sensitive function of the 
stellar parameters that it is justified to treat it as a random quantity. 

In the simplest nonlinear theory of stellar pulsations without resonances, 
sufficiently high amplitude of one mode may stabilize all other linearly unsta- 
ble modes so that they are not present in oscillations. This simple saturation 
mechanism leads to monomode pulsations which are typical for RR Lyrae stars. 

3. Resonant coupling and instability of monomode pulsations 

One of the conditions for the resonant coupling between different modes may 
be written as J2i m i = 0, where rrn are azimuthal numbers of the interacting 
modes. For the most important case of i = 1 it means two possibilities: 

• coupling of the radial mode with anm = nonradial mode, 

• coupling of the radial mode with an m = 1, m = — 1 pair. 

It turns out that the two cases are dynamically equivalent, i.e., the amplitude 
equations describing the mode amplitudes are in both cases the same with the 
same coefficients and in the latter case A± t i = A\-\. 

Radial monomode pulsations become resonantly unstable when 



(Van Hoolst et al., 1998), where k^n is a nonlinear growth rate of the nonradial 
mode, R is a resonant coupling coefficient, Aq^q is a radial mode amplitude, and 
the detuning parameter is given by 



4. Finite amplitude development of resonant instability 

In order to simplify the problem we made some assumptions. The most impor- 
tant are the use of adiabatic resonant coupling coefficients and the neglect of 
nonresonant nonlinear frequency changes. The former assumption was used by 
Van Hoolst et al. (1998) and Dziembowski & Cassisi (1999). The latter assump- 
tion seems to be justified because the nonlinear frequency change has the same 
sign for both modes and it does not change the frequency difference significantly. 

We first studied double-mode stationary solutions of amplitude equations 
together with their stability with respect to small perturbations. The parameters 
of the equations were calculated for the range of stellar models by Dziembowski 
& Cassisi (1999). A typical solution is presented in Fig. 1. All interesting pairs 
of modes in all models give qualitatively the same solutions. 

The a-branch solution exists if the detuning parameter is smaller than the 
critical value given by Eq. (1) (about 0.6 in the case presented in Fig. 1) and it is 
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Figure 1. Double- mode stationary solutions for the radial fundamen- 
tal mode (thick line) and the nearest £ = 1, m = nonradial mode (thin 
line) for a chosen RR Lyrae model. The dimensionless detuning param- 
eter Act is defined as Aw/ko,o- The problem is symmetric with respect 
to the change of the sign of the detuning parameter and the presented 
range of this parameter corresponds to the maximum value given by 
the density of the nonradial modes spectrum. Dotted lines indicate 
unstable solutions. For Act smaller than approximately 0.6 there is a 
pair of solutions denoted by a and b. 



always stable. If the detuning parameter is higher than the previously mentioned 
critical value the monomode solution (not shown in Fig. 1) is stable, according 
to Eq. (1). This means that there always exists at least one stable fixed-point 
solution, either monomode (large Acj) or double- mode (small Aw). 

Direct time integration of the amplitude equations shows that the time- 
dependent solutions always converge to one of the stable fixed-point solutions. 

Multimode fixed-point solutions of resonant amplitude equations are char- 
acterized by the phase-lock phenomenon (for the general discussion see Buchler 
et al., 1997). In our cases it means: 

• for interaction with an m = nonradial mode, the two frequencies are 
equal, 

• for interaction with an m = ±1 pair, the three frequencies are equidistant. 
5. Observational consequences 

When resonant interaction leads to excitation of an m = nonradial mode, an 
observer would see only one frequency due to the phase-lock phenomenon. The 
only effect of the presence of the nonradial mode is the aspect-dependence and 
thus the randomness of the observed pulsation amplitude. 
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The excitation of the m = ±1 pair leads to the observed triplet-type mod- 
ulation, where the side peaks in the spectrum are rotationally split m_ = ±1 
modes. A simple estimate yields the splitting to be equal to £1/2, where Q is the 
mean rotation rate weighted with the Brunt- Vaisala frequency which is strongly 
peaked in the deepest part of the radiative envelope. Thus, in this model the 
Blazhko period is the measure of the rotation in this region. 

6. Conclusions and future work 

The resonant model presented in this paper predicts many observed features of 
RR Lyrae stars: random occurrence and ranges of the modulations as well as 
lower occurence rate of Blazhko stars among RRc than among RRab stars. 

However, this model has some problems, too. The most important one 
is that it does not predict the strong asymmetry of the side-peak amplitudes 
(Kovacs, these proceedings). Moreover, it does not explain features which have 
been recently observed: larger frequency splittings among RRc than among 
RRab stars (Moskalik & Poretti, these proceedings) and the change of the 
Blazhko-period (Smith et al., these proceedings). 

To understand the nature of the long-term modulations of RR Lyrae stars 
more efforts should be made both theoretically and observationally. The most 
promising are attempts to find the nonradial mode characteristics in line-profile 
variations of RR Lyrae itself (Kolenberg et al., these proceedings). Precise spec- 
troscopic observations of other Blazhko stars would also be of great importance. 

Our model also needs some improvements. The most important one is to 
take into account previously neglected coupling coefficients and more interacting 
modes, in particular the whole i = 1 triplet. 



This paper summarizes the results given by Nowakowski & Dziembowski 
(2001) 

Acknowledgments. We would like to thank the organizers for the financial 
support. We are also greatful to Wojtek Dziembowski and Geza Kovacs for their 
useful suggestions given during the preparation of our talk and this paper. 

References 

Buchler, J.R., Goupil, M.-J., & Hansen, C.J. 1997, A&A, 321, 159 

Dziembowski, W.A., Cassisi, S. 1999, Acta Astron., 49, 371 

Kovacs, G. 2000, in Nonlinear Studies of Stellar Pulsations, eds. M. Takeuti and 

D. Sasselov, (Kluwer Academic Publishers), 257, 61 
Nowakowski, R.M. & Dziembowski, W.A. 2001, Acta Astron., 51, 5 
Shibahashi, H. 2000, in ASP Conf. Ser., Vol. 203, The Impact of Large Scale 

Surveys on Pulsating Star Research, eds. L. Szabados & D. Kurtz, (San 

Francisco ASP), 299 
Van Hoolst, T., Dziembowski, W.A., & Kawaler, S.D. 1998, MNRAS, 297, 536. 



